Long-Range Navigation on Complex Networks using Levy Random Walks 
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We introduce a strategy of navigation in undirected networks, including regular, random, and 
complex networks, that is inspired by Levy random walks, generalizing previous navigation rules. 
We obtained exact expressions for the stationary probability distribution, the occupation probability, 
the mean first passage time, and the average time to reach a node on the network. We found that 
the long-range navigation using the Levy random walk strategy, compared with the normal random 
walk strategy, is more efficient at reducing the time to cover the network. The dynamical effect of 
using the Levy walk strategy is to transform a large-world network into a small world. Our exact 
results provide a general framework that connects two important fields: Levy navigation strategies 
and dynamics on complex networks. 

PACS numbers: 89.75.Hc, 05.40.Fb, 02.50.-r. 05.60.Cd 
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I. INTRODUCTION 



Networks are ubiquitous in almost every aspect of the 
human endeavor, explaining the recent burst of work in 
this area [3-13l ■ Besides the topology of the different kinds 
of networks, the dynamical processes that take place on 
them is of utmost importance. In particular, random 
walks are the natural framework to study diffusion, trans- 
port, navigation, and search processes in networks, with 
applications in a variety of systems, such as the propa- 
gation of epidemics and traffic flow [H, 0], animal 
and human mobility T3~ ll[ j an d the dynamics on social 
networks ESS. 



The problem of random walks on networks has been 
addressed before, using a strategy of navigation that con- 
sidered the motion to nearest neighbors [2l|. Here we 
introduce a generalization of this navigation rule by con- 
sidering that the transition probability is not restricted 
to nearest neighbors, allowing transitions that follow a 
power law as a function of the distance (integer num- 
ber of steps) between nodes. This generalized navigation 
rule was inspired by the study of Levy flights where the 
random displacements I obey asymptotically a power-law 
probability distribution of the form P(l) ~ l~ a 22j. 



Random walks on networks are related to the problem 
of searching since one strategy of search is precisely to 
navigate the network, starting from a source node, us- 
ing a random walk until finding a target node [5]. The 
problem of searching and foraging has received consider- 
able attention recently [23 - 25 1 . In particular, it has been 
shown that, under some general circumstances, Levy 
flights provide a better strategy to search or navigate, 
compared with a strategy based on Brownian motion 
24. \26 L For instance, in the problem of foraging by ani- 



mals [24|, ]2]\, |28[ , Levy strategies are rather common, as 



be more efficient than the normal random walk strategy 
(NRW) to cover the network. 

This generalized navigation strategy can consider some 
common situations encountered in real networks. For in- 
stance, in social networks one can take advantage of the 
full or partial knowledge of the network beyond our first 
acquaintances. Currently, using social network sites, one 
can identify the friends of your acquaintances (second- 
nearest neighbors) or the friends of the friends of your 
acquaintances (third-nearest neighbors) and so on, to 
search for a job, an expert opinion, etc. In this way, 
one can contact a second- or third-nearest neighbor of 
a friend directly, without the intervention of the friend. 
This situation corresponds to a long-range navigation on 
a network: a social network in this case. 
Thus, in this paper we are connecting two important 
fields: Levy navigation strategics and dynamics on com- 
plex networks. 

We study the dynamics on an undirected network by 
means of a master equation that we solve exactly with- 
out any approximation. We obtain exact expressions for 
the stationary distribution, the random walk centrality 
iH, the mean first passage time (MFPT) [23], and the 
average time to reach any node on the network. We use a 
formalism in terms of eigenvalues and eigenvectors of the 
transition probability matrix associated with the process. 



II. LEVY RANDOM WALKS ON NETWORKS 

We consider an undirected connected network with N 
nodes i = 1, ...,iV, described by an adjacency matrix A 



with elements A. 
i, j and A4, 



13 ~ 

A; 



well as in the problem of human mobility and behavior 
In a similar fashion, we are proposing that 
our Levy random walk navigation strategy (LRW) can 



Aji = 1 if there is a link between 
= otherwise. We consider the 
case where An — to avoid loops on the network. The 
degree of the node i is given by ki — Ylr—i An. Another 
matrix associated with the network is the distance matrix 
D, with elements dij that denote the integer number of 
steps of the shortest path connecting node i to node j. 



2 



For undirected networks D is a symmetric N x N matrix. 
The average distance (d) scales as a power of N for large- 
world networks, whereas (d) scales as logiV for small- 
world networks 

We start with the discrete time master eq uation that 
describes a random walker on a network 



35] 



A' 



(i) 



where Pij (t) is the occupation probability to find the ran- 
dom walker in j at time t starting from i at t = 0. The 
quantity is the transition probability to move from 
i to j in the network. In the case where i ^ j, it is given 
by ' 



(2) 



where djj are the elements of the distance matrix D 
and Wi->i = 0. This transition probability represents a 
dynamical process where the random walker can visit 
not only nearest neighbors but nodes farther away in 
the network. However, the farther away they are, the 
less probable the event of hopping to that node is. The 
power-law decay of this probability is controlled by an 
exponent a, which is a parameter in our model and 
varies in the interval < a < oo. 

There are two important limiting cases: In Eq. ([2]l. when 
a — > oo we obtain w^j — Aij/ki, which corresponds 
to the normal random walk on networks, previously 
studied by other authors [21], describing transitions 
only to nearest neighbors with equal probability, that 
is, inversely proportional to the degree of the node. On 
the other hand, when a = 0, the dynamics allows the 
possibility of hopping with equal probability to any node 
on the network; in this limit, w^j = (1 — Sij)/(N — 1), 
where 5y denotes the Kronecker delta. For < a < oo, 
the random walker could hop not only to nearest neigh- 
bors but to second-, third- and m- nearest neighbors with 
a transition probability that decays as a power law. 
It is worth noticing that there are navigation strategies 
where some links are added to a d-dimensional lattice 
with a probability proportional to r~ Q , where r is a 
metric distance, generating in this way a small-world 
network j36|. However, in our case, the network is 
not modified, and the possibility of long-range steps is 
determined dynamically by in Eq. ([2]). 

Let us now solve this problem, starting with the sta- 
tionary distribution. By iteration of Eq. (p}, P%j[t) takes 
the form 



E 

h,—,h- 



u lt-l->3 



(3) 



Defining the quantity D 



(a) 



from Eq. ([2]), we 



obta 



Using this relation in Eq. ([3]), 



the detailed balance condition is obtained: 



(0 



(4) 



In order to obtain the stationary distribution, we need 
to take the limit of the occupation probability when the 
time tends to infinity. In this limiting case, the in- 
formation of the initial condition is lost, and the sta- 
tionary distribution gives the asymptotic probability to 
be in a particular node. Thus, for the stationary dis- 



tribution P°° 
r)( a ) poo _ r)( a ) pc 

i 3 j i 



limt^oo Pij(t) [37J, Eq. © implies 



therefore 



D 



(a) 



(5) 



Thus, we have obtained in Eq. (J5j) the exact expression 
of the stationary distribution P°° , which is proportional 
to the quantity given by the sum of the inverse 

of the distances, weighted by a, to node i. Equation 
([5]) generalizes previous results, and it is valid for any 
undirected network. Again, we have two limiting cases: 
When a -> oo, Eq. ^ gives P°° — ^\ ■ , which is 

precisely the result obtained for normal random walks 
yHHI; for a — 0, we obtain the expected result P°° = -h. 

The quantity can also be written as 



JV-l 



(a) = y J_ („) = IH_ >H_ 



(6) 



n=l 



(n) 

where k\ is the number of n-nearest neighbors of node 
i. This equation gives a more clear interpretation of the 
quantity D^: For node «, it is the sum of the first- 
nearest neighbors plus the second-nearest neighbors di- 
vided by 2 Q and so forth. We refer to this quantity as 
the long-range degree. 

We study now the random walk centrality and the 
MFPT. The occupation probability JV, (i) in Eq. (QJ can 
be expressed as 2l|, |35 1 



Pij(t) =S t0 S i: 



Y,PrAt-t')Frj{t') 



(7) 



where Fij(t) is the first-passage probability starting in 
node i and finding node j for the first time after t 
steps. Using the discrete Laplace transform /(s) = 



Ei=o e ~ st /K) in E q- ©. we have 

= \Pi3^)-Sij] /PjM 



(8) 



In finite networks, the MFPT, defined as the mean num- 
ber of steps taken to reach node j for the first time, start- 
ing from node i (3||, is given by (TV,-) = J2tLotFij(t) = 



3 



-F!j(0). Using the moments Pj™ 3 = £^ i™ { P y(*) 
Pj°}, the expansion in series of Pij(s) is 

1 



Pv(s) = P? 



(1 



Introducing this result in Eq. 
MFPT the expression: 

\ 6 J / poc 



(|5J|, we obtain for the 



(0) 



R 



(0) 



(10) 



In Eq. (fTUfl there are three different terms: The mean 



first return time (Ta) 



poo 3 



the quantity Tj 



which is independent of the initial node, and -pis-. The 

j 

time Ti is interpreted as the average time needed to reach 
node i from a randomly chosen initial node of the net- 
work. The quantity Cj, = r^ 1 is the random walk cen- 
trality introduced in 21] and gives the average speed to 
reach the node i performing a random walk. In our case, 
this quantity corresponds to a generalization of the ran- 
dom walk centrality based on the long-range degree given 
in Eq. ©. 

Using Eq. ([TO)) , we obtain 



N 



N 



(ii) 



m— 1 



where the time (T) is the average of the MFPT over 
the stationary distribution P°° and, as shown here, is a 
constant independent of i. In the context of Markovian 
processes it is known as Kemeny's constant and is related 
to a global MFPT j^. 

In order to calculate r, and (T^-) we need to find Pij(t). 
We start with the matrical form of the master equation 
P(t) = P(0)W*. The transition probability matrix W 
is a stochastic matrix with elements Wi-, 



w 



and 

P(t) is the probability vector at time t. Using Dirac's 
notation, 

P ij (t) = (i\W t \j), (12) 

where {|m-}}m=i represents the canonical base of WL N . 
Due to the existence of the detailed balance condition, 
the matrix W can be diagonalized, and its spectrum 
has real eigenvalues [j^. We find a solution of ((T|) in 
terms of the right eigenvectors of the stochastic matrix 
W that satisfy W |<^) = A* |&) for i = 1,..,N. The 
set of eigenvalues is ordered in the form Ai = 1 and 



1 > A 2 > .. > A 



N 



> 



-1. Using the right eigenvectors, 



we define the matrix Z with elements Zij — (i\4>j). 
matrix Z is invertible, and a new set of vectors (<; 
obtained by means of = (4>i\j). Thus 



JY 



The 

>i\ is 



(13) 



Using the diagonal matrix A = diag(Ai, .., Ajv), we obtain 
W = ZAZ" 1 . Therefore, Eq. (JT2J takes the form 



P y (t) = (i|ZA*Z- 1 |j)=^Af(# i )(^li) • (14) 



i=i 



From Eq. ([13]). P9° = (j\(f>i) ((f>i\j) , where the result 

(i\4>i) = const was used. Using the definition of Jw? , we 
have 



N i 

;=2 ! 



(15) 



Therefore, the time Tj is given by 

and for i ^ j in Eq. JTOJ, the MFPT (I^j) is 



i (j|^)<^|j>-(#^)(^b'> 



;=2 



A, 



Wi>(&Ij) 



(17) 



Finally, using ([TT|) and ([T5|) . we obtain 

N N AT 

(T) = EE T3X ^ ; ' TO ) = E t^x • (18) 

m=l !=2 Z=2 ' 

Therefore, we have obtained exact expressions for the 
occupation probability (t) , the stationary distribution 
P°°, the time r i; the MFPT (T y ), and the time (T) in 
terms of the spectrum and the left and right eigenvectors 
of W. Notice that the time (T) is a constant that can be 
calculated using only the spectrum of eigenvalues of W. 



In order to analyze the navigation of a Levy random 
walker on the network, we introduced another global time 
that is simply the average of the quantity n over all the 
nodes on the network, defined as r = jj J2m=i T ™- This 
global time r gives the average number of steps needed 
to reach any node on the network independently of the 
initial condition. In the particular, yet important, case of 
regular networks (such as lattices and complete graphs) 
where the stationary probability is a constant given by 
Pf° = l/N, we obtain that both global times (T) and r 
are the same, r = (T). For example, when a = 0, the 
matrix W has eigenvalues Ai = 1 and Xj — — N _. for 

2,3,..., N, and Eq. JTSJ gives r - i-LzHl 



N 



III. ANALYSIS OF THE RESULTS 

In what follows we use the exact results obtained by 
our matrix formalism, given in Eqs. (|14l) - (|18[) . to cal- 
culate the corresponding quantities. It is important to 
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FIG. 1. (Color online) Ti vs hi for a scale-free network with 
N = 5000 nodes, using three values of the exponent a. The 
inset depicts P°° vs ki , and the solid line indicates the limiting 
case a = 0. 



stress that the results in Figs. [T] and [2] are not numerical 
but exact calculations using the eigenvalues and eigen- 
vectors of the matrix W. In Fig. [T]we show the quantity 
Ti, which represents the average time needed to reach 
node i from any node in the network. We use in Fig. [T]a 
scale-free network of the Barabasi- Albert (BA) type, in 
which each node has a degree that follows asymptotically 
a power-law distribution p(k) ~ k~° |40j. We show three 
cases with different values for the exponent a in Eq. ([2]) , 
that correspond to three different navigation rules. The 
cases a = 1,2 correspond to a LRW, and the limiting 
case a — >• oo corresponds to a NRW 2l| . In the inset we 
show the stationary distribution P°° for different values 
of degree fcj. From Fig. Q]we can conclude that (1) ac- 
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5000 



FIG. 2. (Color online) The time r vs a for different networks 
with N = 5000. (top) The results for a ID lattice, a 2D 
square lattice (50 x 100), and a random tree (network without 
loops); we used periodic boundary conditions in both lattices, 
(bottom) The results for a SF network and an ER network at 
the percolation threshold. 



cording to the inset, using LRW, it is more probable to 
reach the nodes with a small degree (which are the ma- 
jority of nodes on the network), compared with NRW, 
and (2) for nodes with a degree lower than a critical de- 
gree, which is the vast majority, LRW can diminish the 
time to cover most of the network, compared with NRW. 
On the other hand, the time to reach those nodes with 
a high degree is lower using NRW than LRW since LRW 
covers more uniformly the entire network, despite the de- 
gree of the nodes. Therefore, if the purpose is to cover 
the vast majority of the network, it is better to use a 
LRW strategy, but if the goal is to reach the few highly 
connected nodes, it is better to use a NRW. 

In Fig. [21 we show r vs a for five different kinds of net- 
works using our exact results. In regular networks [one- 
dimensional (ID) and two-dimensional (2D) lattices] t is 
calculated using Eq. (fT8)l . and for the other three net- 
works t is obtained by averaging the quantities given by 
(fT6| . The top panel of Fig. [2] shows that for large- world 
networks, such as lattices and trees (networks without 
loops [l|), the LRW strategy navigates the network more 
efficiently. That is, for smaller values of a the average 
number of steps tends to the value (N — l) 2 /N, whereas 
for larger values of a (corresponding to NRW) the num- 
ber of steps can be one or two orders of magnitude larger. 
In the bottom panel we notice that even for small-world 
networks, such as scale-free (SF) and Erdos-Renyi (ER) 
[j| networks, the number of steps is larger for NRW than 
for LRW. Thus, the LRW strategy reduces the global 
time r compared with the NRW strategy, transforming 
dynamically a large- world network into a small world. 

It is worth noticing that the walker can in princi- 
ple visit any node of the network in a single step, but 
with a decreasing probability depending on the exponent 
a. The dynamics involve small, intermediate and large 
steps, drawn from a probability distribution, given by 
Eq. |2]), that decays as a power law. Thus, the dynamics 
is more complex than in a complete graph, where every 
node is connected to every other node; however, as dis- 
cussed before, in the limiting case when a tends to zero, 
the transition probability tends to a constant value given 



by Wi 



(1 — 5ij)/(N — 1), where N is the number of 



nodes on the network. In this limiting case, we arrive 
at a situation similar to a complete graph where each 
node on the network is connected to every other node. 
We show in Fig. [21 that the case a = is optimum in 
the sense that the average global time is minimum; how- 
ever, even for values of a ~ 1, the difference between the 
LRW strategy and the dynamics for a — is negligible. 
This implies that the overall effect of navigating using 
Levy random walks can be as effective as navigating on 
a complete graph. 

These results are confirmed in Fig. [3] where we show 
now Monte Carlo simulations for the percentage of dif- 
ferent visited sites N v in the network, as a function of 
time, for different values of the exponent a. Small values 
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FIG. 3. (Color online) Monte Carlo simulation of the number 
of distinct visited sites N v vs time in a network with N = 
5000. (a) A ID lattice, (b) a random tree, and (c) a scale-free 
network. Each curve is obtained from the average of 1000 
realizations of the random walker. 



(a = 1, 3) correspond to LRW and large values (a = 10) 
to NRW. In Fig. OJJa), we depict the results for a ID lat- 
tice, showing that the number of visited nodes for NRW 
grows diffusively, whereas for LRW grows ballistically. In 
Fig. Efb), we show the results for trees and we see that 
once again the LRW explores the network faster than the 
NRW. Finally, in Fig. [3tc) we show that even for scale- 
free networks there is some advantage in exploring the 
network using the LRW strategy than the NRW strat- 
egy. Even though the difference is smaller, we notice 
that we can cover the network faster using LRW than 
using NRW. 



CONCLUSIONS 

In summary, we have introduced a new strategy of nav- 
igation in general undirected networks, including com- 
plex networks, inspired by Levy random walks, that gen- 
eralized previous navigation rules. We obtained exact 
expressions, using a matrix formalism, for the station- 
ary probability distribution, the occupation probability 
distribution, the mean first passage time and the aver- 
age time to reach a node in any undirected finite net- 
work. We found that using the Levy random walk navi- 
gation strategy we cover more efficiently the network in 
comparison with the normal random walk strategy. For 
small- world networks we obtained that the average time 
to reach any node on the network is smaller than the 
time required using the normal random walk strategy. 



For large- world networks, this difference can be one or 
two orders of magnitude. Additionally, we found that for 
large- world networks, like trees or lattices, the Levy nav- 
igation strategy can induce dynamically the small-world 
effect, thus transforming a large- world network into a 
small world. Finally, our exact results provide a general 
formalism that connects two important fields: Levy nav- 
igation strategies and dynamics on complex networks. 
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CONACYT Mexico. 



[1] M. E. J. Newman, Networks: An Introduction (Oxford 
University Press, Oxford, England, 2010). 
A.-L. Barabasi, |Nature Phys. 8, 14 (2012)1 
S. Boccaletti, V. Latora, Y. Moreno, M. Chavez, and 
D.-U. Hwang, |Phys. Rep. 424, 175 (2006) 



[2] 
[3] 



[8. 
[9 

[io; 
in 
[12: 
[13; 
[14; 
[is; 
[ie; 
[17; 
[is; 

[19 
[20 
[21 

[22 

[23 



and 



A. Arenas, A. Di'az-Guilera, J. Kurths, Y. Moreno, 

C. Zhou, |Phys. Rep. 469, 93 (2008)| 

A. Barrat, M. Barthelemy, and A. Vespignani, Dynami- 
cal Processes on Complex Networks (Cambridge Univer- 
sity Press, Cambridge, 2008). 
A. Vespignani, |Nature Phys. 8, 32 (2012) | 
G. Ramos-Fernandez^ J. L. Mateos, O. Miramontes, 
G. Cocho, H. Larralde, and B. Ayala-Orozco, 
|Behav. Ecol. Sociobiol. 55, 223 (2004)] 

D. Boyer, G. Ramos-Fernandez, O. Miramontes, J. L. 
Mateos, G. Cocho, H. Larralde, H. Ramos, and F. Rojas, 
|Proc. R. Soc. B 273, 1743 (2006)| 

D. Boyer, M. C. Crofoot, and P. D. Walsh, 
| J. R. Soc. Interface 9, 842 (2012) | 

D. Brockmann^ H Hufnagel, and T. Geisel, 
|Nature (London) 439, 462 (2006)] 

M. C. Gonzalez, C. A. Hidalgo, and A.-L. Barabasi, 



Nature (London) 453, 779 (2008) 



M. Bogufia, D. Krioukov, and K. C. daffy, 
|Nature Phys. 5, 74 (2009)] 

G Song, T. Koren, FT Wang, and A.-L. Barabasi, 
|Nature Phys. 6, 818 (20To)] 

V. Belik, T. Geisel, and D. Brockmann, 

|Phys. Rev. X 1,011001 (2011)| 

F. Simini, M. C. Gonzalez, A. Maritan, and A.-L. 
Barabasi, |Nature (London) 484, 96 (2012) | 
D. Lazer, A. Pentland, L. Adamic, S. Aral, et al., 
|Science 323, 721 (2009) | 

J. L. Iribarren and E. Moro, 

|Phys. Rev. Lett. 103 038702 (2009)] 

M. Szell, R. Sinatra, G. Petri, S. Thurner, and V. Latora, 

|Sci. Rep. 2, 457 (2012)1 

S. Aral and D. Walker, |Science 337, 337 (2012)| 
P. Holme and J. Saramaki7 |Phys. Rep. 519, 97 (2012)] 
J. D. Noh and H. Rieger, 

|Phys. Rev. Lett. 92, 118701 (2004) 
R. Metzler and J. Klafter, 

|J. Phys. A: Math. Gen. 37, R161 (2004) 



M. F. Shlesinger. p. Phys. A: Math. Theor. 42, 434001 (2009) 
E. P. Raposo, S. V. Buldyrev, M. G. E. 
da Luz, G. M. Viswanathan, and H. E. Stan- 

R. 



ley, ibid. 42, 434003 (2009) 



Metzler, T. Koren, 



6 



[24] 

[25] 

[26] 

[27] 
[28] 

[29] 

[30] 

[31] 
[32] 

[33] 



B. van den Broek, G. J. L. Wuite, and M. A. Lomholt, 
|ibid 42, 434005 (2009)| G. Osh anin, K. Lindenber g, 
H. S. Wio, and S. Burlatsky, |ibid. 42, 434008 (2009) | 
G. M. Viswanathan, M. G. E. da Luz, E. P. Raposo, 
and H. E. Stanley, The Physics of Foraging (Cambridge 
University Press, New York, 2011). 

O. Benichou, C Loverdo, M. Moreau, and R. Voituricz, 
|Rev. Mod. Phys. 83, 81 (2011)| 

M. A. Lomholt, K. Tal, R. Metzler, a nd K. Joseph, 

|Proc. Natl. Acad. Sci. 105 11055 (2008)] 

D. W. Sims et a?., |Nature (Lo ndon) 45 1, 1098 ( 2008) 
M. de Jager, F. J. Weissing, P. M. J. Herman, B. A . 
Nolet, and J. van de Koppel, |Science 332, 1551 (20lT)| 

C. Brown, L. Liebovitch, and R. Glendon, 
1291 



Human Ecology 35, 
I. Rhee, M. Shin, S. Hong, K. Lee, S. J. Kim, and 
S. Chong, IEEE /ACM Trans. Netw. 19, 630 (2011)] 
N. ScafettaT TChaos 21, 43106 (20l"Tj]~ 
F. Radicchi, A. Baronchelli, and L. A. N. Amaral, 
|PLoS ONE 7, e29910 (2012)] 

F. Radicchi and A. Baronchelli, 

|Phys. Rev. E 85, 061121 (2012)| 



[34] S. Redner, A Guide to First-Passage Processes (Cam- 
bridge University Press, New York, 2001). 

[35] B. D. Hughes, Random Walks and Random Environ- 
ments: Vol. 1: Random Walks (Oxford University Press, 
USA, 1996). 

[36] J. M. Kleinberg, |Nature 406, 845 (2000)] S. Carmi, 
S. Carter, J. Sun, and D. ben Avraham, 

G. Li, S. D. S. 
E. Stanley, and 



Phys. Rev. Lett. 102, 238702 (2009) 
Havlin, 



Reis, A. A. Moreira 
J. S. Andrade, 
Y. Hu, Y. Wang 



H. 



Phys. Rev. Lett. 104, 018701 (2010) 
Li, 



D. 



S. Havlin, and Z. Di, 



Phys. Rev. Lett. 106, 108701 (2011)| 
[37] IfAjv ^ -1, Pi" = lim^oo Pij{t); \i\ N = -1 the Marko- 

vian process is cyclic and P°° is interpreted as a temporal 

average of Pij(t) when t — > oo. 
[38] J. G. Kemeny and J. L. Snell, Finite Markov chains (Van- 

Nostrand, Princeton, 1960); Z. Zhang, A. Julaiti, B. Hou, 

H. Zhang, and G. Chen, |Eur. Phys. J. B 84, 691 (2011)] 

V. Tejedor, O. Benichou, and R. Voituriez, 

|Phys. Rev. E 80, 065104 (2009")] 
[39] N. G. van Kampen, Stochastic Processes in Physics and 

Chemistry (North Holland, T he Netherlands, 1992). 
[40] A.-L. Barabasi and R. Albert, [ Science 286, 509 (1999)| 



